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for obtaining a physical representation of Fourier transform of an 
arbitrary function by means of a torsional device, 

A comparison of such a device is made with two other machines 
known to be in use for obtaining such a representation. It is shown 
that the capabilities of a torsional type transformer are sufficient to 
justify further investigation including the design and construction of a 
practical machine. The basic advantages of a torsional type trans- 
former are simplicity, low cost, and inherent accuracy. 

An investigation of errors arising from the approximation of 
Fourier transforms of certain functions by a step-wise summation 
instead of a continuous integral process is included. The results are 
presented for special cases of interest. 
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OBJECT 



The object of this thesis is to investigate a method for 
obtaining the Fourier transform of an arbitrary function by means 
of a torsional device. 
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CHAPTER 1 



INTRODUCTION 



The development of complex physical systems has many phases and 
sub-phases before the final product is achieved. A brief description of 
these phases, and the inter-relationships between them, may be found 
in the introduction to Volume I of Instrument Engineering (1). Much of 
the information necessary for a successful end product must be generated 
during the actual development. This information may be broad in 
scope and variety, especially in new product developments. It may range 
from physical dimensions and dimensional tolerances to performance 
characteristics of components or the overall system. This thesis de- 
scribes a device which can be useful in certain processes often used to 
obtain performance characteristics experimentally. 

Testing of equipment to obtain actual performance characteristics 
is necessary for several reasons. Although performance characteristics 
may be predicted from theory, experimental verification is needed. 

In other situations testing may prove to be the only means of obtaining 
the necessary design information. 

To this end, transient tests have been performed for many years(2, 3). 
During the early history of closed chain systems correct interpretation 
of the transient responses was difficult. To aid in this interpretation 
frequency analysis methods were borrowed from the practices of 



13 



communication engineering (4, 5, 6, 7). These frequency analysis 
methods include, as the experimental technique, the frequency response 
method. This procedure consists of recording the response of the 
system under test to a sinusoidally varying input. Comparison is then 
made of the steady- state amplitude and phase relationships existing 
between the input and the output. The tests are repeated for the range 
of frequencies desired to obtain the response characteristics as a 
function of frequency. 

It is often inconvenient, and sometimes impossible, to perform 
steady-state sinusoidal frequency tests. An important parameter in 
aircraft dynamic stability calculations is the density altitude at which 
the tests are to be conducted. Hence the frequency method as applied 
to an aircraft in a gliding maneuver would introduce at least one extra 
variable into the test data. Certain components of expendable missiles, 
whose designed working life is only a few minutes, could not be tested 
satisfactorily by tests which would take hours to cover the desired 
frequency range. Other limitations of the sinusoidal frequence tests 
include the possibility of saturation of the equipment in the upper fre- 
quency range. 

It has been shown by Bromberg (8), Seamans, Blasingame and 
Clementson (9) that the frequency response characteristics can also 
be obtained by considering the transient response of the system, partic- 
ularly for a pulse function input. Floyd (10) and Bromberg (8) have 
described the inverse process of obtaining the transient response for 
any forcing function from the known steady- state sinusoidal frequency 
response. 

The relating function, as defined by Draper, McKay and Lees (1) 
concerns the interrelationship between the steady-state frequency 
response, impulse function transient response, and the linear differential 
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equation which is associated with the physical system. The relating 
function concept states that these mathematical notions represent dif- 
ferent ways to describe the basic functional dependence between two 
variables. Special cases of the relating function have been previously 
identified under various names by other authors (4, 6, 11, 12). 

The relating function may be given in several forms. One such 
form is in terms of the relationship existing in the time domain between 
the input variation and the output variation. The relating function may 
also be expressed in terms of the relationship existing between these two 
variables in the frequency domain. When both the input and output vari- 
ations possess zero initial conditions the relating function may be defined 
as the ratio of a response function to a forcing function. The response 
and forcing functions may be obtained as the Fourier transforms of the 
time -dependent response and input functions, respectively. 

Fourier transforms have been derived for many common functions 
and are used extensively in design work. When the input and output 
variations are not available as analytic expressions one of two methods 
must be employed. Either the function must be approximated by a 
function or series of functions for which the transforms are known, or 
an exact method must be found for obtaining the relationship between 
the time and frequency dependent functions. 

The first method presents several possibilities. The function may 
be approximated by a series of rectangular pulse functions of varying 
amplitudes. These pulses are of equal duration and occur at successive 
increments of time. The pulse duration must equal the time increment 
between initiation of successive pulses. Since the transform of the 
rectangular pulse function has been tabulated, the approximate trans- 
form of the original function may be found by a summation process. 
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Series of trapezodial functions have also been used to approximate the 
original function. These methods and their use are discussed in detail 
in Instrument Engineering , Chapter 25 (1). 

Many sets of orthogonal functions may be used to describe variations 
found in nature (13). Of these the most commonly used are the Fourier 
series, which may be expressed as a summation of sines and cosines 
of harmonic frequencies. Fourier series may be- used to determine the 
steady- state sinusoidal relating function if the input and output variations 
contain only frequencies which are integer multiples of a first-order 
frequency. Since the relating function is unique the Fourier trans- 
formation is not necessary if Fourier series representation is used. 
Although in theory the use of a Fourier series to determine the relating 
function would result in an exact solution, in practice only an approxi- 
mation is obtained, since the series contains an infinite number of terms. 
Extensive treatments of the mathematical properties of Fourier series 
are given in many texts (14, 15, 16 ). 

When the steady- state sinusoidal relating function is known and it 
is desired to obtain the response of the system to an input known only 
as a function of time the inverse transformation must be used. This 
procedure is defined as the inverse Fourier transform. With the usual 
assumptions made in engineering practice the approximation methods 
for obtaining the inverse transform closely parallel those for obtaining 
the direct transform. 

The procedure for transforming many functions is laborious and 
time-consuming unless some means of mechanizing the process can 
be found. Machines for this purpose have been constructed (8, 17) and 
have been successfully used. They do, however, possess several basic 
shortcomings. "Devices which approximate the function to be transformed 
by a series of transformable functions have proved to possess limited 
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accuracy at higher frequency ranges. Harmonic analyzers are not 
capable of producing continuous frequency functions. Many, if not all, 
of both types of machines are expensive. 

An investigation has been made into a proposed method of obtaining 
an approximation of the Fourier transform or inverse Fourier transform. 
The device for performing this belongs to the first class of machines 
mentioned, i.e. , those machines which perform the transformation of 
a series of functions. The term "torsional Fourier transformer" has 
been applied to the device, since the principal element is a torsional 
spring. Advantages of the torsional transformer appear to be continuous 
frequency information with greater accuracies and at lower costs than 
machines now in use. 
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CHAPTER 2 



THE FOURIER TRANSFORM 



2 . 1 Definition of a Fourier transform 

The Fourier transform, denoted by the symbol [FT], is defined 
in its most general form by the equation 

[ft] F(x) =fF(>oe~ J **dx =G(Z) 

in which j represents the complex quantity |/- 1 . In practice the equation 
usually indicates a transformation from the time domain to the fre- 
quency domain. Hence the relationship becomes 



[FT] Fit) = J*F<t)<Z~ Jwt dt = GM- 



It should be noted that certain restrictions exist on the original 

function, F(t). If J |F(t) | dt is finite^ and F(t) is continuous, then 

-co 

G(W) exists. The function is usually limited in engineering practice 
to include only functions of the form 



Fit) =0 t<o 
Fit) X o 0<t<T 
F(t)=0 t >T 
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The transformation in this special case has the form: 

[f r] Fit) =// F(t)e ” f dt 

These re strictictions, although limited mathematically, include most 
functions encountered in engineering practice and simplify the problem 
of mechanizing the equation. 

2, 2 The Inverse Transform 

The inverse transformation, from the frequency domain to the 
time domain, is defined as the inverse Fourier transform, represented 
by the symbol [ft] *, and defined as 

JFTY'gM -- Ffr /cWe'H Av = F(t) 

- oo 

Additional conditions which will simplify the equation for the inverse 
transform may be found in Chapter 4. 

This thesis is a description of the design and analysis of a 
device capable of obtaining either the Fourier transform or the inverse 
Fourier transform of a function which meets the conditions outlined 
above . 

No attempt is made in this thesis to define the general properties 
of a Fourier transformation, nor the techniques involved in the solution 
of the transform equation for analytic functions. These features have 
been extensively covered in the literature (12, 14, 15). It should be 
pointed out, however, for those who may be more familiar with the 
LaPlace transformation, that the Fourier forms are special cases of 
the more general LaPlace integrals, in which the parameter jw replaces 
the complex parameter s=c<-f-jw. 
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The Relating Function 



2 . 3 



Potential applications of Fourier transforms for engineering 
purposes have been increased in recent years with the introduction of 
pulse testing of complete physical systems. The relating function 
concepts of Draper and his associates (1) utilize the Fourier transform 
as one of the major tools of analysis and design of linear systems* 

The general form of the relating function is defined by the 
governing equation of the system 



A special case of the general relating function can be defined as the 
quotient of the Fourier transform of the output quantity divided by the 
Fourier transform of the input quantity. This special case will be the 
only case considered in this discussion and will be referred to by the 
more general term "relating function". From the definition of the 
Fourier transform it is seen that this form of relating function is a 
function of frequency; or, as it is often used, a function of a non- 
dimensional frequency ratio. 



linear differential equation this steady-state sinusoidal relating function is 
an operational form of the differential equation relating the input or 
independent variable and the output or dependent variable. The relating 
function, however obtained, is unique and can be used to obtain the 
steady- state response of the system to any sinusoidally varying input. 
Since the theory of Fourier series enables any arbitrary input to be 
expressed in terms of component sinusoidal frequencies, the relating 
function then becomes the only function necessary to express the 
complete behavior of a linear system. 




For any physical system which may be associated with a 
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Applications of the Relating Function 



2. 4 

The implications of this concept have been exploited by Draper, 
McKay and Lees fl), Seamans and others (9, 18), particularly in the 
field of performance function analysis. This particular phase of instru- 
ment engineering deals with the problem of determining the relating 
function of an operating component or system for which the response 
has been determined for a known input. If desired it is even possible 
to obtain from the relating function, determined only as a frequency 
dependent plot, the linear integr o-differential equation which most 
nearly represents the component or system. 

Particular application of pulse methods of testing, together 
with the associated performance function analysis, have been made to 
the field of aircraft stability and control. (9, 18) A pulse input to the 
aircraft is applied in the form of a control surface movement. Recordings 
of both the input pulse and the resulting response of the aircraft are 
used to determine the frequency characteristics. The technique has 
been extended to include the determination of the actual stability de- 
rivatives of the aircraft. Results have compared favorably with those 
achieved as a result of more expensive and time-consuming sinusoidal 
excitation tests. In addition, tests can be performed by the pulse 
testing method under conditions of turn, bank, climb, or glide where 
the sinusoidal excitation technique cannot be used. 

2. 5 Other Applications of Fourier Tran sforms 

Although results have not been published it is known that pulse 
testing methods have been used in evaluating the performance of heat 
exchangers, with results comparable to those obtained by using a 
sinusoidal input. The application of the pulse method to any rate process 
involving linear systems should be straightforward. 
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Fourier transforms of arbitrary functions are also employed in 
seismic research (17). The input pulse consists of a dynamite charge; 
the output is the recorded trace of the earth's vibrations. Since the 
vibrations in which interest is centered lie within a definite frequency 
spectrum Fourier transformation of the time dependent trace is used to 
provide the desired frequency information. The variable area Fourier 
analyzer described in Chapter 3 was developed for this purpose. 
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CHAPTER 3 



MACHINE METHODS OF OBTAINING FOURIER TRANSFORMS 



Computers capable of yielding the transform of an arbitrary 
function fall into two broad classifications, analogue and digital. Although 
the instrumentation differs widely within the analogue group, the governing 
equations for the machines follow one of two basic forms. An analogue 
transformer of each type will be discussed in some detail. Discussion 
of a torsional type transformer will be taken up in Chapter 4. 

3. 1 The Transient Analyzer 

The first of the analogue devices is the Transient Analyzer 
developed in the Instrumentation Laboratory at M. I. T. (8). The analyzer 
is an electro-mechanical machine incorporating 24 resolvers as the 
principal elements. These resolvers are used as variable transformers. 
The rotors of these resolvers are connected to each other by gearing 
arranged in such a manner that the gear ratio between successive 
pairs of rotors is increased by unity. Thus the effective gear ratio 
between the first and last rotors is 24:1. 

This analyzer can be used to obtain the approximate Fourier 
transform of an arbitrary function in the time domain. In general the 
function must satisfy the conditions as outlined in Chapter 2, i.e. , it 
must be a pulse of finite time duration and pulse strength. The period 
of the function is first divided into 24 equal time increments and the 
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amplitude of the function at some average value in each increment is 
determined. This results in the function effectively being replaced by 
a series of step functions of varying amplitude. 

A voltage proportional to the amplitude of the step function 
for the first increment of time is applied to the rotor coils of the first 
resolver. Similar voltages, proportional to the amplitudes of succeeding 
step functions, are applied to the respective rotors. Initially the rotors 
are aligned such that within each resolver the rotor coils are parallel 
to one of the two stator coils. The summation of the voltages from this 
set of stator coils would then be proportional to the integral of the step- 
function approximation. The second set of stator coils, which are 
arranged in space quadrature to the first set, would have no voltage 
output. 

Rotation of the geared rotors by a fixed amount would cause 
any one of the first set of stator coils to have a voltage output equal to 
the original value multiplied by the cosine of the angle through which 
the particular rotor had turned. Similarly, the other stator coil for 
that resolver would have a voltage output equal to the original value 
multiplied by the sine of the angle of rotor rotation. 

That the process results in an approximate Fourier transform 

can be demonstrated. Define F(iAt) as the step-function approximation 

to F(t). Then, since rotation through an angle can be expressed mathe - 

jwAt 

matically as multiplication by the kernel e , one approximation to 

the Fourier transform is given by 

% [F(cat)z- jWiA %t 
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Application of DeMoivre's theorem gives the expression 
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lFT/ Q pp roK F(t) = X F(iAt)[co5 (wt a£)~ j sin (wlaFja t 

Thus the output of the transient analyzer consists of voltages proportional 
to the two summations. 

The error arising from the approximation can be divided into 
two parts. The major error is due to the step-wise, rather than 
continuous, rotation of the function. Errors due to the step function 
approximation of the original function appear to be minor in comparison. 
Further discussion of the magnitudes of the errors will be found in 
Appendix B. For the present it is sufficient to note that this error 
will decrease as the number of time increments increases. Increasing 
this number, however, requires the addition of more resolvers; a 
rather expensive matter. Further complications arise in designing the 
required gear trains for the larger number of resolvers. 

An example of results obtained by the use of the Transient 
Analyzer is shown in Figs. 1 and 2. 

3. 2 A Variable Area Type Fourier Analyzer 

Another approach to the mechanization of the Fourier transform 
is reported by McDonal (17), as applied to problems in seismic research. 
This device was developed for use directly with variable area films 
obtained from a special seismic recorder, although its use with any 
type of variable area film is possible. 

The function, F(t), is converted from a recorder area to a 
voltage by means of a light beam passing through the moving film onto 
a photocell. The amplified photocell voltage is fed into two potentiometers 
whose shafts are turned through angles proportional to cos wt and 
sin wt, respectively. This harmonic rotational motion is obtained by 
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Fig. 1 Phase angle function of Fourier transform of displaced cosine pulse funct ion. 
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Fig. 2 Magnitude function of Fourier transform of displaced cosine pulse function. 



an arrangement of Scotch yokes and pulleys, driven through gear boxes 
by the same constant speed motor which drives the variable area film. 

Thus at any time the outputs of the two potentiometers are 
quantities proportional to the products F(t) sin wt and F(t) cos wt, 
respectively. Integration is accomplished by means of two velocity 
type servos utilizing servo motors and tachometric feedback. The 
outputs, expressed as shaft rotations, are proportional to the two 
components of the Fourier transform for a particular value of the 
variable w. For every additional frequency for which it is desired to 
evaluate the Fourier transform it is necessary to change the gear ratios 
between the driving motor and the Scotch yokes. As constructed the 
device is capable of a frequency range of 150 to 1. 

A fundamental difference exists between the Transient Analyzer 
and the variable area type transformer. The governing equation for 
the latter is 

lFT\F(t) = ff,6)[cos (wt) - j sin (w't’jjdt 

This is seen to be the exact equation for the Fourier transformation, 
and hence no theoretical error would exist. However , the Transient 
Analyzer can be made to yield the approximate frequency function 
directly and continuously, while the variable area type transformer 
would produce only one point of the frequency function for each gear 
ratio used. 

It should be noted that the variable area device has the property 
of producing the Fourier series function corresponding to the original 
time -dependent function. If successive integral values of gear ratios 
are used, successive transformations of a function will result in a 
harmonic analysis of that function. 
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Published reports indicate that a high degree of accuracy can 
be obtained with a device of this type. The disadvantages appear to be 
high initial cost and the labor required to cover a complete frequency 
spectrum. 

3.3 Digital Computers 

Since the physical details of digital computers vary widely no 
attempt will be made to describe an actual machine. The equations 
governing the process of obtaining a Fourier transform are the same 
for all digital computers 

n 

[F F] approK F(t) = X F(c'At)[cos (w cat) -j sin(wiefys,t 

Since this is also the governing equation for the Transient 
Analyzer, the theoretical erros are also the same*. However, a 
decreasing of the time interval, At, would not require an increase in 
equipment for the digital computer. Hence, any desired degree of 
accuracy can be attained. 
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CHAPTER 4 



TORSIONAL REPRESENTATION OF A FOURIER TRANSFORM 



4. 1 Torsional Analogy 

The basis for the torsional analogue of a Fourier transform 

is the analogy between the physical angle of twist of a torsion spring 
—jwt 

and the kernel, e , which appears in the transform equation. Mathe- 
matically, the kernel indicates that a unit length is to be rotated through 
a variable angle in the complex plane. It is necessary to relate the 
physical variables of the spring to the mathematical variables, w and t. 

To relate these variables, consider the effect of twisting a 
long shaft of uniform cross section, holding one end stationary. If 
restraint is not placed on the shaft to retain the original length the 
stresses in the shaft have been shown (19) to be pure shear. Furthermore, 
since the shaft is of uniform cross section, the stress distribution along 
the axis of the shaft will be uniform. This indicates that the angle of 
twist per unit length will be constant. 

Thus the angular displacement of any cross section along the 
shaft is a linear function of two variables the angle through which the 
free end of the shaft has been twisted and the linear distance from the 
cross section to the fixed end. By relating the variable w to the angle 

of twist, and the variable t to the linear distance the analogy of cross 

—jwt 

section angular displacement to the kernel, e , becomes evident. 
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The integrand appearing in the Fourier transform, expression 

- jwt 

is the product of a function, F(t), and the kernel e . Therefore, it 
is necessary to modify the cross section of the shaft in a manner 
analagous to the variation of F(t). This however, must be accomplished 
without changing the linearity of twist per unit length of the shaft. To do 
this exactly would require the addition of material which would impose 
no additional resistance to shear in a plane perpendicular to the axis of 
the shaft, but which possesses rigidity in all planes containing the shaft 
axis. 

This condition can be approximated by adding narrow strips 
of any rigid material mounted on the shaft with their long dimensions 
perpendicular to the shaft axis. Provided that the strips are completely 
separated from each other the condition that there be no restraint in 
shear is approximately true. The degree to which the approximation 
is satisfactory is dependent upon the width of the individual strips. 

Again considering the shaft axis as the axis of the running 
variable, t, the lengths of the strips can be considered as representing 
the function F(t). Figure 3 illustrates such a representation in which 
the function is a displaced cosine pulse. Since F(t) is restricted to a 
real function the function displayed in Fig. 3 must lie in a real plane. 
Therefore, the plane perpendicular to the real plane which contains the 
shaft axis must represent an imaginary plane. This is the familiar 
"contour" representation of conformal transformation theory. 

Rotation of the shaft will result in a presentation similar to 
that of Fig. 4. Projected areas will be seen in both the real and 
imaginary planes. It is the integration of these areas which will result 
in the two components of the desired frequency function. The governing 
equation for the process is 

r? 

[F 7j 0/ r>p rox F(t) = Fftjfcos (wcAt) -j s In (wUtjjA t 
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Fig. 3 Projection of displaced cosine pulse function on real plane (W - 0). 
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Fig. 4 Projection of displaced cosine pulse function on imaginary plane (W /0). 




This approximation differs from that for the transient analyzer 
in that the actual function F(t) is used in the torsional representation 
instead of a step-function approximation. Qualitatively this will reduce 
the error arising from the approximation, although a quantitative analysis 
as given in Appendix B indicates that the major source of theoretical 
error is the same in the two systems. 

4. 2 Inverse Fourier Transforms 

As defined in Chapter 2 the inverse Fourier transform is 
found by evaluating an integral of a function in the frequency domain. 

[FT] f dw = F (F 

In general it can be assumed that G(w) will be a complex function with 
both magnitude and phase angle functions. In engineering practice 
functions most frequently encountered are capable of undergoing further 
simplification. Conditions on the function specify that the magnitude of 
the function G(w) be an even function with respect to positive and negative 
values of w, and that the argument or phase angle be an odd function. 

With the further stipulation that the function does not exist for negative 
values of t, it can be shown (1) that the reduced form of the equation 
becomes 

IrrfcM = /?e.[GMcos(wt)]dw 

where Re is the real part of the function G(w). 

The problem of obtaining the inverse Fourier transform by a 
torsional device is similar to that for the direct transform. The function 
is plotted, either in terms of w or a nondimensional frequency ratio. 
Division of the function into increments of width Aw or is made, 
and the function placed on the torsional shaft. Various values of t are 
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introduced by means of shaft rotation, with integration of projected 
areas being accomplished as for the direct transform. Note that the 
process involves only integration along the real plane, since the 
conditions specified for the function G(w) are equivalent to specifying 
that the inverse transform of G(w) be a real, and not a complex, 
function of t. 

The complete similarity of form which exists between the 
direct and inverse transforms indicates that, by the interchange of 
physical variables, any device capable of performing the direct trans- 
formation is also capable of yielding the inverse transform. Because 
of this similarity an analysis of the inverse transformation process is 
not included in the descriptions of the transient analyzer and the variable 
area type analyzer. 
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CHAPTER 5 



EXPERIMENTAL INVESTIGATION OF A TORSIONAL FOURIER 

TRANSFORMER 



5. 1 Purpose of the Investigation 

The purpose of this investigation was to determine the possi- 
bilities and limitations of obtaining the Fourier transform of an arbitrary 
function by means of a torsional analogue. A description of the equipment 
used to perform this investigation may be found in Appendix A. 

5.2 Selection of Functions 

Functions whose Fourier transforms can be calculated directly 
were chosen to test the performance of the torsional Fourier transformer. 
To carry out the relating function concept functions representing a possi- 
ble input to a physical system and the response of a known system to 
this same input were transformed. 

To this end a displaced cosine pulse was chosen to represent 
the input to a second order system. (Fig. 3). This particular pulse 
shape is particularly adapted to actual testing of physical systems 
since at the beginning of the pulse both the value of the function and its 
derivative are zero. The rectangular pulse function, on the other hand, 
is discontinuous in value at the initial time as well as possessing a 
first derivative of infinite magnitude. As a result it is impossible to 
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build testing equipment for many physical systems which will produce a 
true rectangular pulse function input. The displaced cosine pulse is 
relatively easy to generate and to apply to most physical systems. 

5. 3 Description of Investigation 

A plot of the function to be transformed was made on corrugated 
paper which was then cut into strips of equal width. These strips were 
mounted on a flat spring, one end of which was clamped securely to a 
mounting frame. The other end of the spring was held by a slotted shaft 
which could be rotated with respect to the frame. The entire device 
was placed in front of a lighted opalescent screen. 

Initially the frame was so aligned that the plane of the corrugated 
strips was perpendicular to the line of sight between the frame and a 
camera mounted some distance away. Precautions were taken to mount 
the camera on the perpendicular bisector of the spring axis at a suf- 
ficient distance to minimize the effects of parallax. The first picture 
was taken with the spring untwisted, representing a zero value of fre- 
quency. (Fig. 3). After an angle of twist had been introduced into the 
spring pictures were taken both with the frame in its initial position and 
after being rotated 90°. The process was repeated for additional angles 
of twist. Figure 4 illustrates a typical result for an intermediate 
angle of twist. 

5. 4 Evaluation of Data 

Enlargements of each picture were made to a standard size, 
with a constant enlarger setting. The occluded areas of the enlarged 
pictures represent the desired components of the transformed function. 
Integration of these areas, with due regard to algebraic sign, was ac- 
complished by cutting out the outline of the projected areas and weighing 
the resultant pieces of photographic paper on an analytic balance. 
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The process was readily adaptable to nondimensionalization. 
The reference value of magnitude was taken as the weight of the area 
on the real plane at zero angle of twist. The relationship between the 
angle of twist, Q , and the nondimensional frequency ratio, ^ , for a 
function of period T, is 



O = n W&~t = wT - w T rz + -ZL 

'ref 



since 



w r - 

T«f 

then 

For the results presented in this thesis the reference time 

was taken as the undamped natural period of a second order system. 

The nondimensional frequency ratio, ^ , then becomes the ratio of the 

forcing frequency, w, to the undamped natural frequency, w . It 

n 

should be noted that this choice of reference time has no significance 
for the transform of a function such as a displaced cosine pulse function 
unless the. ratio of the pulse duration to the undamped natural period is 
specified. 



5. 5 Results of the Investigation 

Results of the investigation are presented in Figs. 5^9. The 
magnitude functions are plotted on logarithmic coordinates, the phase 
angle functions on a semi -logarithmic basis. Results of the trans- 
formation of the displaced cosine pulse function are also shown on a 
non-logarithmic plot. (Figs. 1 and 2) 

Included on the figures for comparison are the calculated 
curves of the magnitude and phase angle functions of the Fourier 
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Fig. 5 Logarithmic plot of magnitude function of the Fourier transform 
of a displaced cosine pulse function. 
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a displaced cosine pulse function input. 
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to a displaced cosine pulse function input* 
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transforms. The deviation from the calculated values of the observed 
results are also presented. The transform of the displaced cosine 
pulse function was obtained from a table of transform pairs (1 ), then 
reduced by a process of nondimensionalization to the form 



27T$ 



[F T]fe (1- cos W t)] = 



The response curves of the second order system were obtained 
from analogue computer data (i ). The theoretical Fourier transforms 
of the response functions were obtained by use of the relating function 
concept, employing the relating functions for second order systems to 
be found in Ref. i. 

Mf w = f r 7}yn\RF] s yshtm 

Since both the relating function and the Fourier transform of 
the input function were available in logarithmic form, only an additive 
process was necessary to obtain both he magnitude and phase angle 
response functions. 

Figures 8 and 9 require additional explanation. The entire 
experimental investigation was carried out using the photographic 
method previously described for the integration of areas. However, it 
was conceived that if a practical instrument were to be built, integration 
would probably be dependent upon the use of an integrating photo-electric 
system. As a study of Fig. 4 will show, and as discussed in Appendix C, 
there are gaps between adjacent strips as seen in the projected plane. 
These gaps do not materially affect the envelope of the areas, but would 
presumably affect the values obtained by photo-electric integration. The 
curves of Figs. 8 and 9, therefore, are the results of eliminating these 
gaps from the weighed areas. 
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A derivation of the errors inherent in an approximation of the 
Fourier transform by a step-wise process such as the one used in the 
investigation is given in Appendix B. A partial evaluation of the 
uncertainties present in this investigation may be found in Appendix C. 
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CHAPTER 6 



UNCERTAINTIES IN A TORSIONAL FOURIER TRANSFORMER 



6. 1 Introduction 

A discussion of the uncertainties existing in a representation 
of a Fourier transform by a torsional device must necessarily depend 
upon the details of the actual device used. Since the experimental 
investigation described in this thesis employed certain techniques which, 
would not be used in a practical system the discussion of possible 
sources of uncertainty will not be limited to the experimental methods. 
A more detailed account of the uncertainties considered in the experi- 
mental investigation is given in Appendix C. 

There are two phases to the torsional representation of a 
Fourier transform. The function first undergoes a step-wise rotation. 
Projected areas are then integrated to produce the desired components 
of the frequency function. Uncertainties exist in both the rotation and 
integration phases. 

6 . 2 Rotational Uncertainties 

The use of a torsion spring as the principal element introduces 
several sources of uncertainty. Nonuniformity of cross section or 
lack of homogeneity of the metal will result in variations in the angle 
per unit length. This effect can be minimized by maintenance of 
closer physical tolerances, but can never be completely eliminated. 
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Local differences in the degree of bonding achieved between the 
spring and the strips results in variations in the angular separation 
between adjacent strips. This effect may be decreased by use of a more 
suitable adhesive, or changing the design to allow a larger bonding sur- 
face on each strip. The use of a flat spring with strips fastened to only 
one side caused the strips, when rotated, to become tangent to a circle 
instead of being rotated about the spring axis. This uncertainty can be 
eliminated entirely by the use of a helical spring as the principal element. 
The strips would then be fastened in such a manner that their axes would 
intersect the axis of the spring. 

6 . 3 Uncertainties of Integration 

In the experimental investigation there were two major causes 
of integration errors. As shown in Fig. 4, twisting the spring results 
in a vertical separation occuring between adjacent strips. This is a 
direct result of using a spring with a rectangular cross-section. Sub- 
stitution of a close coiled helical spring for the flat spring would elimi- 
nate most, if not all, of this error. 

The photographic technique used in evaluating the experimental 
data was the source of the other integration error. As explained in 
Appendix A, a weighing method was used to obtain values proportional 
to the integrated areas. Uncertainties arose as a result of nonuniformity 
of paper density. These errors would be eliminated by the use of a 
photo-electric integration system. However, it is believed that the 
photo-electric system would posse ss uncertainties of comparable 
magnitude . 

6 . 4 Summary 

Uncertainties in the torsional transformer representation of a 
Fourier transform occur in terms of angular deviations and integration 
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deviations. By proper design some of the uncertainties present in the 
experimental investigation may be eliminated entirely. All of the un- 
certainties considered may be reduced in designing a more practical 
device . 
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CHAPTER 7 



CONCLUSIONS AND RECOMMENDATIONS 



7. 1 Conclusions 

This thesis was undertaken to investigate the possibilities of 
a torsional Fourier transformer. Experimental investigation has shown 
that a torsional device will produce an approximation to the Fourier 
transform of certain functions. Sources of uncertainty which may be 
present in a torsional representation have been investigated, and 
design criteria established to reduce the level of these uncertainties. 

The two fundamental analogue methods of obtaining a Fourier 
transform approximation have been compared. The basic differences 
have been shown to be 

(1) The method exemplified by a harmonic analyzer will 
produce a transform with no theoretical error, but will 
not provide continuous frequency information. 

(2) The approximation produced by a device such as a torsional 
Fourier transformer contains an inherent error, but is 
capable of being presented as a continuous frequency 
function. 

Inverse transformations have been considered. The ability of 
a torsional transformer to perform the inverse transformation has been 
shown to be subject to the same limitations as exist for the direct trans- 
formation. 



53 



An approximation of the errors inherent in any device which 
obtains the Fourier transform by a step-wise process has been found 
for cases of special interest. The effect of increasing the number of 
steps or increments in the process has been demonstrated. 

7. 2 Re c ommendation s 

Based upon design criteria established in this thesis it is 
recommended that an effort be made to produce a torsional Fourier 
transformer of a practical nature. Specifically, the design should 
incorporate 

(1) a helical spring as the principal element of the torsional 
transformer, 

(2) a photo-electric integration system, 

(3) means for recombining the two component functions to 
obtain a single continuous frequency-dependent function. 

Preliminary studies were conducted to determine a suitable photo- 
electric integration system. Results indicated that a system employing 
two light sources and one photo-cell should be simpler and more accurate 
than other configurations considered. The two light sources, for positive 
and negative area integration, respectively, would be alternately pulsed. 
This system would avoid the necessity for employing balanced photo- 
cell circuits or mechanical shutter arrangements. 

It is suggested that further work be done to determine the in- 
herent errors of the torsional analogue. The possibility exists that 
gain compensation might be applied to reduce the output error of a 
device such as the one suggested above. Such compensation could 
easily be introduced by the same control as used to introduce twist 
into the torsion spring. 
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APPENDIX A 



DESCRIPTION OF THE EXPERIMENTAL TORSIONAL 
FOURIER TRANSFORMER 



The experimental device used to conduct this investigation is shown 
in Fig. 10. There are two separate pieces of equipment, the light box 
and the mounting frame for the torsion spring. Figure 10 shows the 
apparatus in position for an initial photograph, with an arbitrary function 
arranged on the spring. 

The lighting box consisted of a wooden frame, 21 in. square, con- 
taining five 25 watt flourescent light tubes, and fitted with an opalescent 
glass front. The sole function of the box was to provide an evenly lighted 
background for photographic integration. 

Channel aluminum was used in the construction of the mounting frame, 
which measured 18 in. high by 20 in. wide. A steel slotted post was 
fitted in the middle of the bottom channel of the frame, with the slot ar- 
ranged at a 45° angle to the frame. A tapered hole :.n the post at right 
angles to the slot provided a means of securing the end of the steel spring. 
In the middle of the top of the frame a brass slotted shaft was fitted in a 
bushing, with the outer end of the shaft connected by a set screw to a knob. 
The brass shaft was spring loaded to permit motion axially along the flat 
spring which stretched vertically across the middle of the frame. Slots 
in both the shaft and post were cut to a width of 0. 008 in. to prevent move- 
ment of the spring within the slots. 

The heart of the torsional device was the flat spring. Various sizes 
were tried, with the final choice being a strip of spring steel l/4 in. wide 
by 0. 006 in. thick. This appeared to offer a reasonable degree of stiffness 
while permitting sufficient twist to adequately cover the desired range of 
twist angles. 
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Fig. 10 Experimental equipment for investigation of a torsional Fourier transform 



Finding a suitable material from which to prepare the function rep- 
resentation presented another problem. Advertising display paper, con- 
sisting of a corrugated layer bonded to a backing sheet, proved to possess 
sufficient rigidity when cut in thin strips. It was found from experience 
that best results were obtained by cutting a rectangular piece of paper into 
l/4 in. strips which were left joined at one end. The spring was then 
fastened to the unjoined ends of the strips, and the function was cut to the 
desired shape. The use of a standard paper cutter is recommended for 
the initial cutting of the strips; the slightly jagged edgeleft by a razor blade 
or scissors results in interference between adjacent strips. 

The most time-consuming problem encountered was the finding of a 
suitable method of securing the paper to the spring. Although the paper 
strips were light in weight, the bonding surface on each strip was only 
l/l6 sq. in. In addition to the weight effects, the twisting of the spring 
resulted in additional shearing forces being exerted on the bond. Of the 
many adhesives tested only one was found to be satisfactory; a double- 
sided pressure-sensitive tape produced by Minnesota Mining and Manu- 
facturing Co. Best results were obtained by taking photographs immedi- 
ately after fastening the paper to the spring. 

It should be noted that it is necessary to mount the spring with its 
axis vertical, if gravity effects on the uniformity of twist are to be elim- 
inated. 

Photos were taken, as described in Chapter 5, with a 35 mm. camera 
mounted at a distance of eight feet from the spring. The camera was 
mounted on the perpendicular bisector of the spring axis, with the line of 
sight also perpendicular to the real plane of the function. 



57 



APPENDIX B 



ERRORS INVOLVED IN APPROXIMATIONS 
TO A FOURIER TRANSFORM 



The error involved in approximating a Fourier transform by a step- 
wise process such as is used in the torsional Fourier transformer is 
difficult to evaluate. Quantitative errors vary greatly with the degree 
of regularity or "smoothness" of the function transformed. Laning has 
shown (20) that the error calculation inevitably leads to an infinite series 
solution, the form of the series being dependent upon the form of the 
approximation of the original function. 

Therefore, emphasis in this appendix is primarily upon evaluating 
the effect on the magnitude of the error function of increasing the number 
of increments of the running variable Since the error function is of 
more general interest when expressed in non-dimensional terms, it is 
necessary to first define the non-dimensional form of the Fourier trans- 
form. 



This is exactly equivalent to the dimensional form of the Fourier 
transform as defined in Chapter l. 

To evaluate the -error involved in an approximation to the Fourier 
transform of an arbitrary function assume an approximation of the 
form 





where 
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(£) 



which is the general case of the approximation inherent in the torsional 
Fourier transformer. 



Define 



(l ) [FT] FCt) = FT] FCt) ~ n rr o, Fft) 



then 
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substituting in (3) and collecting terms 
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At this point it becomes necessary to make certain assumptions 
concerning the function, F('t). One reasonable assumption is that 



F('t)~ f (i A 't) -h ''T [FUA 't)J 



within the interval iA't < 't < (i + 1 ) A 't . 
Substitution into (6) yields 



(£)(FT]FC6)= I F(LA’t)<z: j21r§LA ' t f%-^' T -i) d r 



C ~ ~oo 



/ i [Fa A n e^Y'rF 7 * ' r -/)d ‘r 



C = -co 



(?) 



Performing the indicated integration 
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Since 



[FUa t)]A 't = F(ch)a 1 - F(Ca ' t) 
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eq. (8) can be expressed as 



(S )[FT] FCt> = 



-jZITQA't 
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In considering only the magnitude of the error function it is seen 
that the first summation can be expressed as the product of a deviation 
factor and a Fourier transform approximation. The second summation 
includes a term involving the difference between the value of F ( 1 1 ) at 
successive points. This latter term, from the conditions imposed upon 
F('t), must change sign during the period of the function. Furthermore, 
for functions which are symmetric, such as the displaced cosine pulse 
function, the effect on the magnitude of the error function disappears. 

Therefore, an approximation to the magnitude of the error function 
becomes 



il){FT\F(‘t) 



(c 



.jznBA't ) 

~ J) -I 



-jzrrgk't 



F T \ rr , x FCt) 



This differs from the usual definition of a deviation function, which 
is defined as the ratio of the error to the true value. This can be ex- 
pressed as 



(D)[FT]FCt) = 



rFTlFCt)-tFrl^ m FCt)_ , (-FWAt l 
[FT] F(’t) 



This deviation function is plotted in Fig. 1 1 for two values of t, 
together with observed deviation curves of an input function transform 
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Fig. 11 Deviation functions of Fourier transform approximations. 






and a response function transform. Several observations should be made, 
however, concerning a direct comparison between the observed and cal- 
culated curves. 

First, for Fig. 11 is taken as the period of the function, hence 

the apparent disagreement between the deviation curve for the response 
function and the deviation curve of Fig. 9, in which T re ^ was taken as the 
undamped natural period of the system. Secondly, the fact that the ob- 
served deviation is less than the theoretical deviation for the input func- 
tion may arise from the method of obtaining the integrated areas for this 
case. The rapid increase in the observed deviation functions near 
$ = 2. 0 is to be expected since the theoretical value of the transforms 
at this point is zero. 

No direct comparison with the transient analyzer results can be 
made, since the transformed function in this case was more nearly 

£■[ F(i+OA 't + F(lA V 
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APPENDIX C 



UNCERTAINTIES IN THE EXPERIMENTAL INVESTIGATION 



The strict definition of the term "uncertainty 11 relates to an inac- 
curacy that is erratic and can only be treated on a statistical bases. In 
this sense at least one of the inaccuracies discussed in this appendix is 
not an uncertainty, since it is amenable to direct calculation. This 
error (item 2 below) is, however, an error in an intermediate physical 
quantity, and will result in an output error which will vary with the func- 
tion being transformed. It is therefore treated as an uncertainty. 

There are three general sources of uncertainty inherent in the ex- 
perimental setup, and one source arising from the method chosen to 
evaluate the data. Since a thorough quantitative investigation of these 
sources was not attempted the listing is not necessarily in order of mag- 
nitude. 

(1) Variation in the angle of rotation existing between adjacent 
strips. 

(2) Uncertainties due to the fact that the function was not rotated 
about the actual axis of the spring 

(3) Vertical separations appearing between adjacent strips as 
projected into the component planes. 

(4) Variations in photographic paper density. 

Items (1) and (2) concern uncertainties effecting an intermediate 
quantity, the angle through which the strips have been rotated. Items (3) 
and (4) are inaccuracies in the integration process which produces the 
output quantity. 

The variation in the angles between adjacent strips was measured for an 
extreme angle of twist. The average angle between strips was 14. 2°, the 
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maximum 16. 5° and the minimum 1.0°. The variation between adjacent 
groups of three angles, however, was much less, the minimum being 
42° and the maximum 45°. This indicates that for a regular function 
such as those transformed most of the individual errors resulting from 
the non-uniform rotation would cancel out. 

The second source of uncertainty listed occurred as a consequence 
of fastening the paper to one side of the spring. As a result the strips 
do not rotate about the axis of the spring, but instead are tangent to a 
small circle. The radius of this circle is equal to the distance between 
the axis of the spring and the axis of a strip. This distance was of the 
order of 0. 02 in. for the particular combination of spring, adhesive tape, 
and paper strips used. 

The gaps appearing between adjacent strips as seen in the projected 
planes are caused primarily by the distortion of the spring. It has been 
demonstrated (19) with a shaft of rectangular cross section that lines 
drawn on the face perpendicular to the axis are distorted when the shaft 
is twisted. Since the strips are fastened along the entire width of the 
spring, this distortion produces a shear force between the spring and the 
strips. The resultant action of the individual strip then becomes depend- 
ent on which part of the distorted face of the spring remains firmly bonded 
to the strip. The magnitude of this uncertainty can be estimated by the 
difference in deviation between the curves of Figs. 7 and 9. Since only 
the envelope of the projected area was used in determining the curves of 
Fig. 7 the effect of the distortion of the spring does not appear. Figure 9, 
however, was the result of eliminating the lighted area which appeared 
between adjacent strips. A direct determination of the uncertainty can 
not be made for two reasons. The transformed functions were not identi- 
cal for the two figures, although they did possess the same general char- 
acteristics. There is also the possibility that part of the difference in 
deviations was due to the errors calculated in Appendix B. 

Uncertainties arising from non-uniform density of photographic 
paper were partially investigated. A critical case, in which the differ- 
ences between the positive and negative areas was small, was chosen for 
this purpose. The results obtained by the weighing method were com- 
pared with those obtained by carefully graphing the photographed areas 
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and integrating by counting squares. Results differed by less than two 
percent of the reference magnitude. 
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